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Abstract 

We present a technique for estimating the number of future Ov/3/3 results using several distinct 
nuclei to optimize the physics reach of upcoming experiments. We use presently available matrix 
element calculations and simulated sets of predicted 0^/3/3 measured rates in multiple isotopes 
to estimate the required precision and number of experiments to discern the underlying physics 
governing the mechanism of the process. Our results indicate that 3 (4) experimental results with 
total uncertainty (statistical, systematic, theoretical) of less than ~ 20% (~ 40%) can elucidate 
the underlying physics. If the theoretical (i.e. matrix element) uncertainty contribution is below 
~ 18%, then 3-4 experimental results of ~ 20% precision (statistical and systematic) are required. 
These uncertainty goals can be taken as guidance for the upcoming theoretical and experimental 
programs. 
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I. INTRODUCTION 



The study of neutrinos in general, and that of double-beta decay ((3(3) in particular, is 
among the most exciting fields on the cutting edge of fundamental physics research as it will 



guide our development of the standard model of partic 



study (3(3 is well motivated in the literature 



ii,y,y,y. 



e physics. In particular, the need to 
5)]. Two-neutrino double-beta decay 



(2is(3(3) (e.g. 76 Ge — > 76 Se + 2e~ + 1v) is a second-order weak decay process, and is the most 
rare nuclear decay allowed in the standard model. 2v(3(3 can occur in even-even nuclei where 
beta decay is forbidden, and has been observed in approximately ten nuclei. Zero-neutrino 
double-beta decay (Qv(3(3) violates lepton number conservation and is therefore forbidden 
in the standard model. 0^(3(3 differs from 2v(3(3 by the absence of antineutrinos in the final 
state (e.g. 76 Ge — > 76 Se + 2e~). Though the observation of Qv(3(3 would imply the existence 
of massive Majorana neutrinos[6 : ] (that is, massive neutrinos that are indistinguishable from 
their own antiparticle), a number of possibilities exist for the underlying physics mediating 
the decay. Figure [U shows the Feynman diagrams for 2u(3(3 and 0^(3(3. 

This article will focus on the use of a suite of experimental results of the 0v(3(3 ground- 
state transitions in different nuclei to distinguish exchange mechanism models or transition 
matrix-element models. This approach is of immediate interest because of the significant 
number of next-generation experiments that are being proposed that might provide these 
data. In particular, motivating the required number of experimental results from different 
nuclei is critical to the overall (3(3 program. 

In Section [IIJ we begin with an overview of the problem of calculating 0i>(3(3 matrix ele- 
ments, discuss the exchange mechanisms considered in this study, and compile a review of 
the current matrix element calculations. Section II III reviews previous work similar to that 
presented here. Section IIVI covers the details of how we performed the model comparison. 
Section IVl estimates the required total uncertainty as a function of the number of experimen- 
tal results. In Section IVIl we summarize and draw some conclusions about this technique's 
impact on the discovery potential of future 0^(3(3 searches. 

II. Ou(3(3 MECHANISM AND MATRIX ELEMENT CALCULATION OVERVIEW 

The Op (3(3 rate can be written: 
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FIG. 1: A hadron- level Feynman diagram of 2uj3/3 (left) and Qv(3(3 (right) for the light neutrino 
exchange mechanism. 



G» v \M 0uV \ 2 (1) 

where T 0u is the decay rate, G 0u is the kinematic phase space factor, \M 0u \ is the 

matrix element corresponding to the 0vf3j3 transition, and rj is a lepton number violating 
parameter (LNVP). G 0u contains the kinematic information about the final state particles, 
and is exactly calculable to the precision of the input parameters (though use of different 
nuclear radius values in the scaling factors of G 0u and \M 0u \ have previously introduced some 
confusion [7]). The LNVP J] contains all of the information about lepton number violation, 
and has a form depending on the 0v(3j3 mechanism. 

The LNVP takes on different forms for different mechanisms. In addition, \Mq v \ also 
depends on the mechanism. In this paper we consider: light Majorana neutrino exchange, 
heavy Majorana neutrino exchange, right-handed currents (RHC), and exchange mechanisms 
that arise from R-Parity violating supersymmetry (RPV SUSY) models. Our model choices 
for consideration in this paper were motivated by the existence of calculated \Mq v \ for many 
isotopes within those models. The heavy-particle models represent a large class of theories 
that involve the exchange of high-mass (>1 TeV) particles. For example, leptoquarks 



very similar |Mj 
current models 



_. have 

to RPV SUSyQ. Left-right symmetric models can lead to right-handed 
lOj or heavy neutrino exchange models 11]. Scalar bilinears[12] might also 



mediate the decay but explicit matrix elements have not been calculated yet. For SUSY and 
left-right symmetric models, effective field theory|l3l| has been used to determine the form 
of the effective hadronic operators from the symmetries of the GV/3/3-decay operators in a 
given theory. This last paper makes clear the close connection between all the heavy-particle 
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exchange models. 



A. Matrix Element Calculations 

The nuclear matrix element in equation [I] contains the information about the nuclear 
structure of the reaction, and is influenced by the exchange mechanism considered for the 
transition. Two formalisms have been used to calculate |M ^|: the shell model and the 
quasiparticle random phase approximation (QRPA). It is important to note the relative 



strengths and weaknesses of the QRPA 



16] 17] and shell model approaches [18} . QRPA 



is 



capable of including an effectively unlimited number of single-particle states, but a limited 
number of configurations. The opposite is true of the shell model. In this light, it is 
encouraging that the two methods arrive at similar answers. Most recent calculations of 
\M 0u \ have been done in the QRPA, renormalized QRPA (RQRPA) or proton-neutron QRPA 
(pnQRPA) frameworks because the nuclei of interest for double-beta decay are sufficiently 
large as to make the shell model calculations numerically complex. Encouraging work, 
however continues in the shell model, and the confidence in this technique will only improve 
as computing speeds and numerical methods become more advanced. Overviews of the 
matrix element calculations are available in the literature [19I. I20I . Recent works by Rodin et 



al 



2l| and Suhonen 2^ 

provide a written debate concerning the implementation of QRPA. 
This debate has helped elucidate why the collection of previous QRPA calculations have 
resulted in numerous differing results. In particular, the free parameters in the theory can 
be normalized differently. Additionally, not all calculations have included the same input 
physics. 

B. Ovf3(3 Mechanisms and the Lepton Number Violating Parameter 

In this section we define the LNVP for each of the models. 



1. Light Neutrino Exchange 

If the neutrino is a Majorana particle, it can be exchanged between two neutrons and 
mediate If the neutrino is light, the LNVP t]l u has the form: 
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£ viz 



m e 



(2) 



where the index k spans the light neutrino states, U e k is the (e, k) element of the neutrino 
mixing matrix, is a phase, is a neutrino mass eigenstates eigenvalue, and m e is the 



electron mass. Detailed discussions o 
element can be found 



in [l6 
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ight neutrino exchange Ov/3/3 and its associated matrix 



21. 



22 



23 



24j . Further description of the interplay 



between the LNVP and the neutrino mixing parameters is in 



2. Heavy Neutrino Exchange 



Heavy Majorana neutrinos can also contribute to 0upf3. The LNVP (rjuv) however, has 
a different form: 



VH, = J2U^. (3) 

Here, k runs over the heavy neutrino states, is the mass eigenvalue of the heavy neutrino 
(mfc >> lGeV), and m p is the mass of the proton. It is important to note that for light 
neutrino exchange, rjL U is proportional to the light neutrino mass. In contrast, r]Hu is inversely 
proportional to the heavy neutrino mass. Heavy neutrino exchange is discussed at greater 



length in Ref. [24j] and references therein. 



3. Right-Handed Currents 

In left-right symmetric models[10(, right-handed currents can result in interactions leading 
to OvPfi. Two parameters are used to described this phenomenon, rj denotes the magnitude 
of a possible right-handed leptonic current coupling to a left-handed hadronic current in the 
weak interaction Hamiltonian. A denotes the magnitude of coupling between a right-handed 
leptonic and a right-handed hadronic current. For more information about A and rj see 
reference 1251. 
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FIG. 2: The two-nucleon or short-range contribution (left), and the long-range In (middle), and 
In (right) exchange mechanisms in the heavy-particle exchange framework. 



4- R Parity Violating Supersymmetry 



Many supersymmetric (SUSY) R-parity violating extensions to the standard model have 
lepton-number violatin g in teractions between quarks and leptons that can also contribute to 



26 



27 



28 



29 



30 



31 



32 



331 ] . The particles exchanged are heavy and the physics takes 



place at short distances. Hence, the heavy degrees of freedom can be integrated out leaving 
hadron-lepton vertices. The two-nucleon (2N) or short-range contribution to Qvf3f3 tends to 
be suppressed because of the large mass of the exchanged particle. Alternatively, the heavy 
particle can convert into a virtual n, where the comparatively low-mass n meson can mediate 
the Qvf3j3 transition between more distant nucleons. If one of the initial quarks gets placed 
into the virtual n, we refer to it as the In mode. Similarly, if both are placed into a virtual 
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271 ] . Generally, 



7r, we have the 2-7T mode. These two modes are discussed in detail in Refs 
many of the heavy particle exchange models have similar \Mq u \ with similar short-range and 
long r ang e components. The direct exchange and its relative suppression are discussed in 



Refs. 
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26. 



27 



28] . Figure [2] shows the Feynman diagrams for these 3 modes. 



For the SUSY 0v(3f3 mechanisms, the product of M and t]susy take the form 26]: 



M VsusY = VT M 2 - N + ( VP s - V r)Mf + ^( V t + ^VpsK^M 1 * + M 2w ) (4) 

Here, rj? and rjps are the effective lepton number violating parameters that normalize the 
tensor and pseudoscalar currents in the lepton number violating part of the effective La- 
grangian. M^ N and are the matrix elements for the two-nucleon direct exchange mode, 
and M l7r and M 2n are the matrix elements for the n exchange mode. 
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C. Matrix Element Tabulation 

In Table HI we tabulate the phase space integrals in units of 10 15 y. G 0u includes a scaling 
factor (r ), which is traditionally factored out and combined with a similar factor from the 
matrix element to conveniently make the matrix element unitless. However a problem arises 
if authors use matrix elements and phase space integrals from different sources and are not 
careful to make certain that they are scaled by the same value of r j3]. In Table [J the 
numerical value of r$ used is tabulated for each reference. 



Reference 


r (fm) 


48 Ca 


76 Ge 


82 Se 


96 Zr 


100 Mo 


116 Cd 


128 Te 


130 Xe 


136 Xe 


150 Nd 


[24] 


1.1 




7.93 


35.2 


73.6 


57.3 


62.3 


2.21 


55.4 


59.1 


269 


m 


1.2 


75.8 


7.60 


33.5 


69.7 


54.5 


58.9 


2.17 


52.8 


56.3 


249 



TABLE I: Phase space factors and the assumed nuclear radius scaling factors. The units are 
10~ 15 /year. 



Table [TT] contains matrix elements for light neutrino exchange calculated in either the 
QRPA or shell model framework. There are literally dozens of other matrix element cal- 
culations that are not tabulated here. As discussed in Refs. {2, 21], many of the previous 
calculations are outdated. We chose the QRPA values for light neutrino exchange in the 
Table, because we consider them to be the state of the art. We include values for the QRPA 
from Ref. 



24 



25| for comparison only, since we use values for heavy neutrino exchange 



and right-handed currents from those references, respectively. For the Shell Model, we used 
the published values from Ref. [18| except for the cases of 100 Mo and 116 Cd, where we have 
taken the liberty of using the values from Ref. |35|] because those values weren't included in 
Ref. [is)]. Since the values in Ref. js^] have not yet been published, one must consider their 
use preliminary and only indicative of the anticipated final results. 

In Table IHIt we tabulate the matrix elements arising from heavy neutrino and SUSY 
particle exchange. The choice of reference for the heavy neutrino exchange mode was made 
due to its estimates of \M 0v \ for numerous isotopes. For the SUSY mode, it is the only 
comprehensive set of calculations done after the realization of the importance of the Itt and 
2ir modes. We assumed that the l-zr and 2ir terms in Eqn. H] dominate, and the values in 
the table are taken from Ref. 26j] as: 
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M SUSY = (-M^ + M^). (5) 



Reference 


Method 


r (fm) 


48 Ca 


76 Ge 


82 Se 


96 Zr 


100 Mo 


116 Cd 


128 Xe 


130 Te 


136 Xe 


150 Nd 


[24] 


pnQRPA 


1.1 




2.8 


2.64 


1.49 


3.21 


2.05 


2.17 


1.8 


0.66 


3.33 




RQRPA a 


1.1 




2.4 


2.12 


0.31 


1.16 


1.43 


1.6 


1.47 


0.98 


2.05 


[17] 


pnQRPA & 


1.2 




3.33 


3.44 


3.55 


2.97 


3.75 




3.49 


4.64 




[25] 


pnQRPA c 






14.1 


12.3 




3.6 




15.2 


9.6 


2.0 


28.9 


[18] 


Shell Model 


1.2 


0.720 


1.39 


2.19 








0.86 


1.19 


0.75 




[35] 


Shell Model 


1.2 




1.6 


1.7 


0.4 


0.3 


1.9 


2.0 


1.6 







a gA = 1.25, fix coupling parameters to reproduce correct 2v(3f3 rate 

6 fix coupling parameters to reproduce single (3 rate of intermediate nucleus 

c This paper expresses the matrix element and phase-space factors as a product and doesn't explicity state 
the r value. 



TABLE II: Qv(3f3 nuclear matrix elements for several light v exchange models. 



Reference 


Method 


r (fm) 


48 Ca 


76 Ge 


82 Se 


96 Zr 


100 Mo 


116 Cd 


128 Xe 


iso Te 


136 Xe 


150 Nd 


[24J 


pnQRPA 


1.1 




32.6 


30.0 


14.7 


29.7 


21.5 


26.6 


23.1 


14.1 


35.6 


[26J 


pnRQRPA 


1.2 




-625 


-583 


-429 


-750 


-435 


-691 


-627 


-366 


-1054 



TABLE III: 0v(3(3 nuclear matrix elements for heavy neutrino exchange (upper row) and R-parity 
violating SUSY models. 



Table ITVl provides the matrix elements arising from right-handed currents as calculated by 
Ref. 251 ] . Again we chose this reference because it provided values for a number of isotopes. 
However, it does not provide a value for 116 Cd. Therefore, in the analysis described below, 
these two theories require special treatment. 



III. SUMMARY OF PREVIOUS WORK 



If 0v(3f3 is observed, it will demonstrate that neutrinos are massive Majorana particles 
regardless of the underlying physics of the process The prospect of finding a distinctive ex- 
perimental fingerprint for this underlying physics would then become extremely interesting. 
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Reference 


Method 


48 Ca 


76 Ge 


82 Se 


96 Zr 


100 Mo 


116 Cd 


128 Te 


i30 Xe 


136 Xe 


150 Nd 


[25] 


pnQRPA(ry) 




0.44 


1.54 




3.50 




0.15 


2.25 


0.83 


29.5 


[25] 


pnRQRPA(A) 




0.14 


1.01 




1.05 




0.007 


1.05 


0.20 


26.8 



TABLE IV: Products of nuclear matrix elements and the corresponding phase space factors 

for right-handed-current interaction models. The units for the rj values are 10 _8 /year and for the 
A values are 10 _12 /year. 



Various techniques have been discussed previously. The use of the kinematic distributions 
to discern right-handed currents (RHC) from light neutrino mass processes was described 



in Ref. [101. Ref. 



361 ] discusses the use of the angular distribution as a function of particle 



physics model. In Ref. 



371 ] . the author also proposes examining the ratio the rates of 0vf3(3 



to the 2 + excited state to that for Qisfifl to the ground state in the same isotope as a sig- 



nature for RHC. Ref. 38j discusses using relative rates of the decay to the first excited + 
state and the ground state to discern the light neutrino mode, the heavy neutrino mode, 
and SUSY processes. Ref. [39j discusses the potential of using fi — > e and [i — >■ e7 in com- 



bination with OuPP to prove that the light v mechanism is dominant. Ref. [40j discusses 
the relative contributions to T 0u from light neutrinos and heavy particle exchange. In this 
report, we revisit the approach of comparing rates of the ground state transition in multiple 
nuclei to determine the mechanism. Singlet neutrinos in extra dimensions can lead 

to double beta decay [2]. In this case the Kaluza-Klein spectrum of neutrino masses spans 
values from much less to much greater than the nuclear Fermi momentum. Therefore one 
cannot separate the nuclear and particle physics and the effective neutrino mass depends 
on the nuclear matrix element. The decay rate also depends on unknown parameters such 
as the brane shift parameter and and the extra dimension radius and therefore is highly 
model dependent. Finally, mass varying neutrinos (MaVaNs) might lead to decay rates that 



depend on the density o 
Three recent articles 



the matter in which the process occurred 



41 



42 



15|- 



43 [have considered the question of the required number 



of 0v(3/3 experiments. Reference 



4l| took a very pessimistic view of nuclear theory and 



PP. That paper took all previous calculations of \Mq u \ and considered them as samples 
of a parent distribution of the true value. This led to a very wide distribution of possible 
anticipated results for a given (m ee ) and hence a requirement that an unpractically large 
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number of measurements be performed. Many of the calculations, however, have been found 
to be in error, or have been updated and hence replaced by more recent calculations. (See 
Ref. {21] for a summary discussion.) If one neglects the outlying results, the situation is 
improved. 

Reference 42j concludes that Qufi/3 measurements in 3 nuclei would be an important tool 



in the solution of the nuclear matrix element problem. This reference does not address the 
question of discerning the various mechanisms for 0u(3j3, although it does point out that 
it is an interesting one. It also does not quantify the confidence one would derive from 3 
measurements nor does it quantify the required precision. 

During the final preparations of this manuscript, Ref. 43j became available. This paper 
estimates the ratio of rates for 0uj3P for a number of mechanisms but doesn't estimate the 
number or required measurements or their required precision. The paper does estimate the 
spread in ratios for a number of particle physics models and suggests pairs of isotopes that 
would be most useful for separating certain models. 



IV. MULTIPLE-ISOTOPE COMPARISON 



To compare theoretical calculations and perspective experimental results, we calculate 
the predicted rate for each considered model using Eqn. [TJ The LNVP and \M 0u \ appear as 
a product. Hence, if one obtains a decay rate from a lone nucleus, the two factors cannot 
be individually separated. That is, two different LNVP values can lead to the same decay 
rate. Therefore, we normalized the LNVPs for each model so that all reproduce identical 
decay rates for 76 Ge. The absolute normalization was chosen for (m ee ) = 100 meV with the 
matrix elements in Ref. [sij]; that is, the scale of rj corresponds to Of/3/3 lifetimes just within 
the reach of the current set of proposed experiments. This absolute scale for the LNVP is 
arbitrary for this analysis and is only chosen for definiteness. Our arguments concerning the 
uncertainty requirements and the number experimental results do not depend upon it. 

Our model space will consist of seven different theories enumerated below chosen from 
those described in Section Hi CI Because the RHC model does not have a \M 0u \ for 116 Cd, 
our primary analysis uses only the first 5 of these models. A secondary analysis tests the 
impact of including the RHC models. 

1. QRPA(/3/3): Light v exchange, RQRPA with coupling constants tuned to reproduce 
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the 2v{3(3 rate [21] 



2. QRPA(/3): Light v exchange, pnQRPA with couplings tuned to reproduce the /3-decay 
rates for the intermediate nucleus in the 2v(3(3 reaction [17J 



3. SM: Light v exchange, shell model la. l35| 

4. Heavy v: Heavy v exchange 24] 



5. SUSY: R-parity violating supersymmetry for the dominant l7r and 2ir exchange 
modes [261] 



6. RHC-??: Right-handed current 25| 



7. RHC-A: Right-handed current [25] 



For each of these models, we simulated 10000 potential sets of decay-rate results for a group 
of isotopes with measured decay rates according to a Gaussian distribution whose mean was 
defined by the predicted decay rate and its variance defined as a fraction of that mean. We 
then compared each simulated result group to each of the models to discern the predictive 
power of selecting the correct model. Specifically, each simulated group was compared to all 
models by calculating x 2 - The best model that fit the simulated set was chosen as that with 
the lowest x 2 value. We then tabulated the number of simulated sets for which the best 
choice was also the correct set. We then determined the largest total uncertainty that would 
result in the correct choice being selected 90% and 68% of the time. This analysis was done 
not only for varying uncertainty, but also for a varying number of isotopes included in the 
analysis. (Alternatively, one could tabulate the number of times a wrong model was chosen 
as the correct model and determine how well one can reject a given model.) 

The six (3f3 isotopes used in this analysis are: 76 Ge, 82 Se, 100 Mo, 116 Cd, 130 Te, and 136 Xe. 
These six were chosen because they are proposed for future experiments and have \Mq v \ 
calculations available for a variety of models. 

The analysis presented below is by no means restricted to this set of Qv[3(3 models and 
isotopes. In fact, one would anticipate that the theories will improve and that experimenters 
will invent new techniques that employ other isotopes. Most importantly, the theory is not 
entirely settled. Therefore, although we use the T 0u values quoted in Table \V\ as indicative 
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of the difference between the various models, we recognize that the predictions are likely to 
evolve. It is clear that this work rests on a number of assumptions, including: 

1. We have made the implicit assumption that the present differences in the model pre- 
dictions are indicative of the true differences. 

2. It is conceivable that future calculations will include additional microphysics (e.g nu- 
clear deformation) that will affect different isotopes differently. Our analysis assumes 
that any uncertainties in the models shift each \Mq v \ similarly. If a model has a sys- 
tematic shift in its predictions for all isotopes that alters all \Mq v \ by an equivalent 
fraction, that would not alter this analysis in choosing the correct theory, but would 
result in a shift in the value of the deduced LNVP. 

3. We have assumed that within each set of theoretical values for |M ^|, each calculation 
for the individual isotopes is of equivalent accuracy. 

4. Finally, we have also assumed that one and only one model is correct. That is, we 
have assumed that if two mechanisms are contributing to T 0u (e.g. light neutrino and 
heavy particle exchange), one mechanism dominates and any interference is negligible. 

These assumptions notwithstanding, the analysis can provide useful guidance indicating how 
precise measurements should be and how many measurements are required. 

V. MODEL SEPARATION ANALYSIS 

From the calculated T ^ we determine the mean values and uncertainties for each isotope. 
The explicit values are plotted in Fig. [3] with 10% uncertainties shown to provide a qualitative 
indication of how well the various models can be discerned at this chosen uncertainty level. 
Values are given in Table |V] This table is very useful for identifying which models will be 
difficult to separate. For example, the SUSY and heavy neutrino models give very similar 
decay rates for each isotope to within a somewhat uniform factor of about 2. Hence it will 
be hard to separate these two models from each other. Also the table shows the interesting 
case of the shell model estimate for 100 Mo. This isotope shows a large disagreement (greater 
than a factor of 10) with several of the other models. Hence, this Table indicates that Mo is 
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FIG. 3: r predictions for \Mq v \ uncertainties of 10% for all 7 Qv(3(3 models and all 6 isotopes. 



a key isotope for separating the models. (However note the caveat regarding the reference 
for the shell model calculations discussed above. This value has not yet been published.) 



Model 


LNVP 


Vb Ge 


82 Se 


iUtJ Mo 


iib Cd 


130 Te 


i3b Xe 


QRPA(/3/3) 


1.96 x 10" 


-UY 


1.75 


6.1 


3.0 


4.9 


4.6 


1.2 


QRPA(/3) 


1.41 x 10" 


-07 


1.75 


8.3 


10.0 


17.4 


13.4 


25.3 


SM 


3.38 x 10" 


-07 


1.75 


19.3 


0.6 


25.7 


9.0 


3.8 


Heavy v 


1.44 x 10" 


-08 


1.75 


6.6 


10.5 


6.0 


6.1 


2.4 


SUSY 


7.51 x 10" 


-10 


1.75 


6.7 


18.2 


6.6 


12.3 


4.5 


RHC-7? 


6.30 x 10" 


-10 


1.75 


6.11 


13.88 




8.92 


3.28 


RHC-A 


1.14 x 10 


-7 


1.75 


13.07 


13.59 




13.59 


2.64 



TABLE V: Values of the decay rate in units of 10 27 /year for the isotopes in our analysis for the 
chosen LNVP. The light neutrino mass ((m ee )) was chosen to be 100 meV for the model labeled 
Rodin. The other LNVP's were chosen such that the resultant decay rates for 76 Ge were the same. 
The form of the LNVP was chosen to be unitless. 



To determine the number of isotopes required to obtain sufficient separation, we varied 
the number of isotopes included in the analysis. By design, the analysis does not include 
all isotopes, and therefore, the order in which isotopes were added to the analysis had to be 
chosen. We chose four isotope orderings; by atomic number, by largest spread in predicted 
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T 0u , by an assumed likely order of actual results, and an alternate ordering to examine 
the effect or 100 Mo. The "experimental readieness" (third) ordering and the choice of the 
isotopes is clearly arbitrary and reflects the opinions, not any prescience, on the part of the 
authors. Note: 76 Ge is first in each list is because, as mentioned in Section HVl we normalized 
the LNVPs to reproduce the same rate in that isotope. 

1. Atomic number: 76 Ge, 82 Se, 100 Mo, 116 Cd, 130 Te, 136 Xe 

2. r°- spread: 76 Ge, 136 Xe, 116 Cd, 100 Mo, 82 Se, 130 Te 

3. Likely order of experimental results: 76 Ge, 130 Te, 136 Xe, 100 Mo, 82 Se, 116 Cd 

4. Alternate order to study 100 Mo: 76 Ge, 130 Te, 136 Xe, 116 Cd, 82 Se, 100 Mo 

All of the (3(3 isotopes treated in this article have proposals for upcoming Fiorini-style 



internal source experiments with the exception of 82 Se. The Majorana 44j and GERDA 
collaboratk 
experiment 



45] 



collaborations will use 76 Ge. The MOON collaboration 46j plans to use 100 Mo. The COBRA 

49] 



47| will use 116 Cd. CUORE|48j will use 130 Te and the EXO collaboration 
will use 136 Xe. 82 Se is one of the (3(3 isotopes proposed for use by the SuperNEMO 
collaborationist]] in a tracking apparatus. 

Table [VTl summarizes the results of the simulations. 



VI. CONCLUSIONS 



If we examine the results in Section El we see that the ordering of the isotopes has a 
noticeable impact on the results. Although, this ordering dependence does indicate that cer- 
tain models might be better tested with particular isotopes, given the assumptions regarding 
the nuclear theory upon which these calculations are based, it would be dangerous to put 
too much weight on the specific choice of isotope for an experiment based on this work using 
presently available matrix elements. When considering only 5 models (rows 1-8 in Table 
IVip . the results imply that approximately 4 experimental results are required to discern the 
best physics model underlying the Qv(3(3 process, depending on the confidence level desired. 
For 68% (90%) confidence in choosing the correct model, a total uncertainty (theoretical 
and experimental) of ~40% (~20%) might be acceptable if 4 experimental results are avail- 
able. In this analysis it is the heavy neutrino model that is the limiting factor in setting the 
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TABLE VI: Limits of required total uncertainties corresponding to each theoretical model and 
the number of measurements as isotopes are added to the analysis for 90% and 68% confidence. 
Only the final pair of rows includes all 7 models. The others exclude the RHC models because 
those models do not have a \Mq u \ value for 116 Cd. The precision of the numbers in this table are 
approximately ±1-2%. 



uncertainty levels. If only light neutrino models are considered, then the requirement for 
68% (90%) can be relaxed to >50% (~35%) uncertainty for just 3 experiments. From this 
analysis, it is also clear that only 2 experimental results is not sufficient. 

When all 7 models are used, the required statistical precision becomes more stringent than 
when one uses a smaller number of models, as one would expect. A 68% confidence- level 
result would require 4 measurements in different isotopes with a precision of ~25% or better. 
(See the final pair of rows in Table I VI I ) The general conclusion remains that separation 
is possible if the uncertainties can be made small, but it also indicates that additional 
constraints limiting the possible viable models from non-Oz//3/3 experimental results will help 
discern the underlying mechanism. 

The uncertainty described above is the total uncertainty: statistical, systematic and theo- 
retical. The experimental uncertainty (statistical and systematic) for the proposals discussed 
above can be expected to reach below 20% or so. This leaves a remaining requirement on 
the theoretical uncertainly of ~35% for 68% confidence-level separation of the 5 models. 
Because \Mq u \ appears squared in the formula for the decay rate, the matrix element theo- 
retical uncertainty is required to be half this value or <18%. This uncertainty requirement 
on the Qv/3P matrix elements may seem somewhat daunting given their historical spread, 



15 



but as discussed in Section 1111} there has been tremendous progress in understanding the 
source of this spread. This is especially true in the case of QRPA-type calculations, pro- 
viding hope that 10 — 15% uncertainties are possibile. In fact, the difference between the 
shell model and QRPA \M Qu \ results are already within ~30-40% for some isotopes. Since 
these two approaches are radically different, it seems reasonable to use this difference as an 
estimate of the precision of these calculations. Such an estimate, however, does not replace 
the need for uncertainties for the calculations derived from theory itself and from auxiliary 
measurements [5 1| . In any case, should calculations reach these uncertainty levels, it will 



greatly increase the physics reach of 0v(3/3, and help to formulate the future program of 
searches for physics beyond the standard model. 



53| . The 



There are other motivations that require multiple 0^/3/3 experimental results 
need to prove that the observation is indeed and not an unidentified background is 

not the least among them. Therefore, a general conclusion from this work, and those similar 
to it, indicate that at least 3 (and very likely 4) 0v(3(3 experiments along with significant 
theoretical effort are warranted. The need for and utility of several precision experimental 
results is the critical conclusion from this work. 
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